Are the anti-charmed and bottomed pentaquarks molecular heptaquarks? 
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I study the charmed uuddc resonance D*~p (3100) very recently discovered by the HI collabo- 
ration at Hera. An anticharmed resonance was already predicted, in a recent publication mostly 
dedicated to the S=l resonance © + (1540). To confirm these recent predictions, I apply the same 
standard quark model with a quark-antiquark annihilation constrained by chiral symmetry. I find 
that repulsion excludes the D*~p (3100) as a uuddc s-wave pentaquark. I explore the D*~p (3100) 
as a heptaquark, equivalent to a N — tt — D* linear molecule, with positive parity and total isospin 
I = 0. I find that the N — D repulsion is cancelled by the attraction existing in the N — tt and 
7r — D channels. In our framework this state is harder to bind than the Q + described by a k — tt — N 
borromean bound-state, a lower binding energy is expected in agreement with the HI observation. 
Multiquark molecules N — tt — D, N — n — B* and N — tt — B are also predicted. 



I. INTRODUCTION 



In this paper I study the anti-charmed uuddc reso- 
nance D*~p (3100) (narrow hadron resonance of 3099 
MeV decaying into a D*~p) very recently discovered 
by the HI collaboration at HERA pj. This extends to 
flavour SU(4) the SU(3) anti-decuplet HHH which in- 



cludes the recently discovered 6+ (1540) |5j, ffl, 0, IE 
and H— (1860) [13, E3, Thc + (154O), S" (I860) 
and D*~p(3100) are extremely exciting states, because 
they may be the first exotic hadrons to be discovered, 
with quantum numbers that cannot be interpreted as 
a quark and an anti-quark meson or as a three quark 
baryon. Exotic multiquarks are expected since the early 
works of Jaffe [IS El EE EE, and the SU(3) exotic 
anti-decuplet was first predicted within the chiral soliton 
model 0, IS Hi- Pentaquark structures have also been 
studied in the lattice 0, 0, 0] . However the nature 
of these particles, their isospin, parity |20j and angular 
momentum, are yet to be determined. 

We recently completed a work on the uudds + |2l| 
where we indicate that the + is probably a K — tt — N 
molecule with binding energy of -30 MeV. In that work 
we first compute the masses of all the possible s-wave 
and p-wave uudds pentaquarks, and we verify that these 
pentaquarks are hundreds of MeV too heavy to explain 
the + resonance, except for the 1 = 0, J p = l/2 + state. 
However in this channel we find a purely repulsive exotic 
N — K hard core s-wave interaction. This excludes, in our 
approach, thc + as a bare pentaquark uudds state or as 
a tightly bound s-wave N—K narrow resonance. We then 
add thc tt — N, tt — K and N — K interactions to study 
the + as a borromean three body s-wave boundstate of 
a n, a N and a K |H |H EE] , with P ositive parity ^ 
and total isospin 1 = 0. In that paper, and in a very 
recent work [23], we also address the S= -2, Q = -2 state 
discovered by the NA49 experiment with a mass of 
1.862 GeV, indicating that this is a K — — K molecule 
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with a binding energy of -60 MeV. 

In ref. [2l| we also conclude suggesting the existence 
of similar anti-charmed uuddc and anti-bottomed ex- 
otic uuddb hadrons. The anti-charmed pentaquark was 
widely expected 22], and its properties will certainly con- 
tribute to clarify the nature of the pentaquarks. Thc 
state D*~p (3100) may be similar to thc + , with the an- 
tiquark s replaced by a c. In this case it is natural to con- 
sider replacing the K meson by a D meson or by a D* me- 
son, because the D* is also a narrow state. For instance 
in the new positive parity D s mesons, [24| the K — D 
and the K — D* multiquarks are respectively candidates 
to thc scalar D s (2320) and to the axial D s (2460). Thc 
energy of thc is D*~p (3100) consistent with a D* — tt — N 
linear molecule with an energy of +15 MeV above thresh- 
old. This case differs from the previous ones because here 
there is no negative binding energy. Nevertheless a sys- 
tem which energy is located slightly above threshold is 
still a narrow state, and in this sense the D*~p (3100) 
remains in the same family of thc + and of thc 5 . 

Moreover a natural theoretical motivation exists for 
considering heptaquarks (or pentaquarks) and not just 
p-wave pentaquarks (or baryons) in the exotic multiplct 
of the + . Supose that a given s-wave pentaquark hadron 
H is studied and one concludes that it is unstable. Nev- 
ertheless one may consider that a flavor singlet quark- 
antiquark pair uu + dd or ss is created in the hadron H. 
When the resulting heptaquark H' remains bound, it is 
a state with an opposite parity to the original H , where 
the reversed parity occurs due to the intrinsic parity of 
fermions and anti-fermions. In this sense the new hep- 
taquark H' can be regarded as the chiral partner of H . 
And, because H' is expected to be aproximately stable, it 
is naturally rearranged in a s-wave baryon in two s-wave 
mesons. The mass of the heptaquark H' is expected to 
be slightly lower than the exact sum of these standard 
hadron masses due to the binding energy. This princi- 
ple explains qualitatively the mass of the + [2l| for the 
S [23j and the masses of the non-exotic multiquarks 
of the D s and D* 8 familly H3- 

In this paper I extend the techniques used in our first 
publications to the D*~p (3100) and to the other simi- 
lar narrow resonances with an anti-charm or anti-bottom 
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quark. A standard Quark Model (QM) Hamiltonian is 
assumed, with a confining potential and a hyperfine term. 
Moreover the Hamiltonian includes a quark-antiquark 
annihilation term which is the result of spontaneous chi- 
ral symmetry breaking. I start in this paper by reviewing 
the QM, and the Resonating Group Method (RGM) 
which is adequate to study states where several quarks 
overlap. Using the RGM, I show that the corresponding 
exotic baryon-meson short range s-wave interaction is re- 
pulsive in exotic channels and attractive in the channels 
with quark-antiquark annihilation. The short range re- 
pulsion contradicts the existence of narrow pentaquarks 
with an anti-charm or anti-bottom quark. I proceed 
with the study of the linear molecules or heptaquarks 
D - tt - N, D* - tt - N, B - tt - N and B* - tt - N. In 
particular the total energy of these systems is discussed. 
Finally I conclude interpreting the D*~p (3100) and pre- 
dicting related multiquarks. 



II. FRAMEWORK 

Our Hamiltonian is the standard QM Hamiltonian, 

E T *+E^+E^ d) 



H 



where each quark or antiquark has a kinetic energy Ti 
with a constituent quark mass, and the colour depen- 
dent two-body interaction Vij includes the standard QM 
confining term and a hyperfine term, 
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(2) 



The QM of eq. reproduces the meson and baryon 
spectrum with quark and antiquark bound-states (from 
the heavy quarkonium to the light pion mass). The RGM 
was first applied by Ribeiro 29] to show that in exotic 
N — N scattering, the quark-quark potential together 
with the Pauli repulsion of quarks explains the N — N 
hard core repulsion. Deus and Ribeiro [3(j also showed 
that, in non-exotic channels, the quark-antiquark anni- 
hilation could produce a short core attraction. Recently, 
addressing a tetraquark system with the tt — tt quantum 
numbers, it was shown that the QM also fully complies 
with the chiral symmetry, including the Adler zero and 
the Weinberg theorem [3lL l32j | . Therefore the QM is ad- 
equate to address the anti-decuplet, which was predicted 
in an effective chiral model. 
For the purpose of this paper the details of the poten- 
tials in eq. are unimportant, only its matrix elements 
matter. The hadron spectrum constrains the hyperfine 
potential, 



(V hyp ) ~ - (M A - M N ) ~ M 
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When a light quark is replaced by a heavy quark, say 
a charmed quark, the hyperfine interaction is decreased, 
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FIG. 1: Examples of RGM overlaps are depicted, in (a) the 
norm overlap for the meson-baryon interaction, in (b) a ki- 
netic overlap the meson-meson interaction, in (c) an interac- 
tion overlap the meson-meson interaction, in (d) the annihi- 
lation overlap for the meson-baryon interaction. 



and it must also be replaced by (Vh yPD ) — Mr>* — Dk 
. The quark-antiquark annihilation potential Aq is also 
constrained when the quark model produces spontaneous 
chiral symmetry breaking |33l l34| . The annihilation po- 
tential A is present in the tt Salpeter equation, 



2T+V A 
A 2T + V 



(4) 



where the tt is the only hadron with a large negative en- 
ergy wave- function, <p~ ~ <p + . In eq. Q the annihilation 
potential A cancels most of the kinetic energy and confin- 
ing potential 2T+V. This is the reason why the pion has 
a very small mass. From the hadron spectrum and using 
eq. @ the matrix elements of the annihilation potential 
arc determined, 



(2T + V) S=0 
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N 



Ma 



\{2M N -M A ) , 



(5) 



where this result is correct for the annihilation of u or d 
quarks. 

The RGM |2^| computes the effective multiquark en- 
ergy using the matrix elements of the microscopic quark- 
quark interactions. Any multiquark state can be de- 
composed in combinations of simpler colour singlets, the 
baryons and mesons. The wave functions of quarks are 
arranged in anti-symmetrized overlaps of simple colour 
singlet hadrons. Once the internal energies Ea and Eb 
of the two hadronic clusters are accounted, 



(^q|g£ P (-l) P ^a< 
(Ma\E P (-^) P P\^l 



E a + E b + V al 



(6) 



where J2 P (~ 1) P -P is the anti-symmetrizer, the remain- 
ing energy of the meson-baryon or meson-meson system 
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is computed with the overlap of the inter-cluster micro- 
scopic potentials, 

V^a = (<P B ^A\-(V li + V 15 + 2V 2i + 2V 25 )3P u 

mcs B 

+3A 15 \cf) A <f> B }/(0 B A |1 - W 14 \<p A <p B ) 

V z ,a = {(f>B<f>A\(l + PAB)[-{Vl3 + V23 + Vu + V 2 4) 

xP 13 + A 23 + A u ]\(j) A (p B ) 

/{cf> B cf> A \(l + P AB )(l - Pn)\M B ) , (7) 

where Pjj stands for the exchange of particle i with parti- 
cle j, see Fig. ^ It is clear that quark exchange provides 
the necessary colour octets to match the Gell-Mann ma- 
trices Xi present in the potential Vij . This results in eq. 
© or eq. JHJ) times an algebraic colour x spin x flavour 
factor and a geometric momentum overlap |3J|. 

A good approximation for the wave-functions of the 
ground-state hadrons is the harmonic oscillator wave- 
function, 



for the different channels, are 



0ooo (Pp) =N a x exp 



2a 2 J 



20F 



(8) 

where the inverse hadronic radius a can not be estimated 
by electron-hadron scattering because it is masked by the 
vector meson dominance, a is the only free parameter in 
this framework. In the case of vanishing external mo- 
menta pa and p Bl the momentum integral in eq. J7J) is 
simply N a ~ 2 ■ 

The annihilation potential only occurs in non-exotic 
channels. Then it is clear from eq. @ that the annihila- 
tion potential provides an attractive (negative) overlap. 
The quark-quark(antiquark) potential is dominated by 
the hyperfine interaction of eq. ||3J|, and in s-wave sys- 
tems with low spin this results in a repulsive interaction. 
These results are independent of the details of the quark 
model that one chooses to consider, provided it is chiral 
invariant. Therefore I arrive at the attraction/repulsion 
criterion, 

- whenever the two interacting hadrons have quarks (or 
antiquarks) with a common flavour, the repulsion is in- 
creased by the Pauli principle, 

- when the two interacting hadrons have a quark and an 
antiquark with the same flavour, the attraction is en- 
hanced by the quark- antiquark annihilation. 

In the particular case of one nucleon interacting with 
anti-kaons and with kaons, this implies that the short 
range exotic D - N, D* - N, B - N and B* - N in- 
teractions are repulsive. This shows that the 1=0 s-wave 
pentaquarks uuddc and uuddb are certainly quite unsta- 
ble. Higher isospin or spin systems are certainly more 
unstable in our framework because they have a higher 
repulsion. On the other hand the short range tt — N, 
tt — D, tt — D* , tt — B and tt — B* interactions can be 
attractive. This motivates the study of a linear molecule 
with a N, a tt and a D, or a D + , or a B, or a B*. Quanti- 
tatively |2ll|32,|36|,|33, the effective potentials computed 
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where f are the isospin matrices, normalized with f — 
t(t+ 1). The vanishing momentum case of eq. © is suf- 
ficient to compute the scattering lengths with the Born 
approximation. However the study of binding needs the 
finite momentum case. In the exotic D — N and D* — N 
channels, it can be proved that the geometric result Af a ~ 2 
is then replaced by the separable interaction |0ooo)( < /'oool- 
In the non-exotic fr — N, tt — D and tt — D* channels the 
present state of the art of the RGM does not allow a 
precise determination of the finite momentum overlap. 
Nevertheless I assume for simplicity the same separable 
interaction. This is a reasonable approximation, because 
the overlaps decrease when the relative momentum of 
hadrons A and B increases. Moreover when the hadronic 
potential V AB is in a separable form v\4>\) {<pi\, the energy 
of the bound state and the matrix element of the poten- 
tial are simple to compute. A boundstate coincides with 
a pole in the T matrix at a negative energy, 
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therefore one just has to find the energy that cancels 
1 — go v. This method also allows the computation of the 
matrix element of the potential in the boundstate, 



< V AB >= v 



go 



E 7 'fo,i 2 



When |0i > is the harmonic oscillator state of eq. 
the necessary condition for binding is, 

— 4/x v > a 2 . 



(11) 



(12) 
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III. BINDING FLAVOUR uuddQ MULTIQUARKS 

The simplest pentaquarks are not expected to bind 
due to the attraction/repulsion criterion. For instance 
the D*~p (3100) cannot be the ground-state uuddc pen- 
taquark because the elementary color singlets (uud)-(dc) 
or (udd) — (uc) are repelled, since the elementary color 
singlets share the same flavour u or d. This also implies 
that the D~p and D°n systems are unbound. The uuddc 
pentaquarks with spin, flavour or angular momentum ex- 
citations will have larger masses and also large widths. 
Nevertheless the c pentaquarks are more subtle than the 
s ones, because the D* is quite stable when compared 
with the K*. Therefore one should also consider to ex- 
cite the spin in the Ic cluster, and this amounts to study 
D* — N bound-states. Indeed the effective potential of 
eq. © is attractive in this case. However this state is 
coupled to the D — N case also in eq. (0 . Once the cou- 
pled channel hamiltonian is diagonalized, the energy of 
the D* — N is lifted and the attraction is essentially lost. 
Therefore the simplest way to have attraction, together 
with a low energy and with a narrow width, consists in 
adding at least one quark-antiquark pair to the system. 

Then this amounts to include a pion in the system, 
which can be attracted both by the N baryon and by the 
D or D* meson to produce a,N — tt — DoyN — tt — D* 
linear molecule. For instance the flavour includes combi- 
nations of terms like uud — du — uc where the anti-quark 
u in the pion can be annihilated both by the u present in 
the N and by the u present in the D. According to the 
attraction/repulsion criterion this produces an attractive 
interaction. The quark d present in the nucleon cancels 
only part of the attraction to the pion. Incidently the 
pion- nucleon I = 1/2 attraction is fixed by chiral sym- 
metry, see reference |36j| . 

The proposed system N — it — D and N — it — D* are 
similar to the model for O + (1540) advocated in reference 
[21). in the present case the anti-quark s is replaced by 
a heavy c or b. The increase of the quark mass does 
not affect directly the attraction, where the Q is just a 
spectator. However the size of the wave-functions 1/a is 
affected. For instance in an harmonic oscillator potential 
a is proportional to ^//Z, and the reduced mass u doubles 
when one changes from a light-light meson to a heavy- 
light meson. This amounts to an increase of nearly 20% 
of the a in the D or B meson. Because the a parameter 
is increased only in one of the Jacobi coordinates, the 
average a in tt — D or tt — D* or tt — B or tt — B* is 
only expected to suffer a 10%. This increase of a will de- 
crease effectively the attractive interaction, for instance 
the condition for a 2-body binding in eq. (|12|) shows 
that this is equivalent to decrease the strength of the 
attractive potential by a factor of 20%. Similar results 
are obtained in different models of confinement, say in 
the funnel interaction which is more adequate for heavy 
quarks. In what concerns the repulsive D — N potential, 
it is decreased in the same way. Moreover the strength 
of the hyperfine potential is further decreased because 



(Vh yPD ) << (Vhyp)- For example the strength of the re- 
pulsive D — N potential is expected to decrease by 30 
%. 

I now use an adiabatic Hartree method to study the 
stability of the linear N — tt — D molecule and related 
molecules with a, D* , & B oi & B* . Essentially the wave- 
function of the pion is centered between the nucleon and 
the -D, where the nucleon and the D don't overlap with 
each other. This results in a linear molecule, for simplic- 
ity I use an averaged mass for the nucleon and D and 
for the pion interaction with these quark clusters. I solve 
a Schrodinger equation for the nucleon in the potential 
produced by a pion placed at the origin and by the other 
heavy-light meson placed at a distance —a of the pion. 
The potential of the pion is produced by the D meson at 
the point —a and the nucleon +a. This produces three 
binding energies EjjE^, En, and three wave- functions. 
In the Hartree method the total energy is the sum of these 
energies minus the matrix elements of the potential en- 
ergies, This is easily computed once the two Schrodinger 
equations are solved, with eqs. (TDJ, 111|) and l|12l) . The 
total energy is a function of the distance a, and I mini- 
mize it as a function of a. The same steps are repeated 
for the N - tt - D* , N - tt - B and N - tt - B* sys- 
tems. At the point I am not yet able to bind these linear 
molecule systems, with a negative binding energy. The 
same happened for the N — tt — K when we studied the 
+ . [21| Nevertheless the picture of a K — tt — N with 
a binding energy of 30 MeV is still plausible because the 
medium range interaction remains to be used. Therefore 
binding or near-binding is also plausible in the N—tt — D, 
N - tt - D*, N - tt - B and N - tt - B* systems. 



IV. CONCLUSION AND OUTLOOK 

I find that N - tt - D, N - tt - D* , N - tt - B and 
N — tt — B* nearly bound 1=0 linear s-wave molecules, 
and positive parity are plausible. The absorption of a low 
energy pion and the resulting decay into a low energy p- 
wave ./V — D(D*, B, B*) results in a narrow decay width. 
The N—TT—D* and N—tt—B* can also decay respectively 
into the three body systems N — tt — D and TV — tt — D, 
but this is again narrow since the N — D* and N — B* 
overlaps are suppressed. Moreover N — D*, and N — B* 
pure 1=0 s-wave, negative parity, pentaquarks are not 
completely ruled out, but they couple strongly with the 
decay channels N — D and N — B (repulsive systems), 
and this results in a larger width. 

Comparing with the K — tt — N and assuming that it is 
a bound system, I expect the N —tt— D* and N — tt — B* 
to be the most bound systems of the family studied in 
this paper, with a similar binding energy. Nevertheless 
the N — tt — D* is less bound than the N — tt — K, due 
to a weaker tt — D(D*, B, B*) attraction when compared 
with the ttk attraction. Thus an energy of 14 MeV above 
threshold, corresponding to a uuddc mass of 3.099 GeV 
as observed by the the HI collaboration is plausible. The 
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corresponding mass of the multiquark with flavour uuddb 
is of the order of 6.416 GeV. I what concerns the gound- 
state N — ir — D and N — tt — B the binding energy is 
predicted to be higher because in this case the nucleon- 
heavy- light meson repulsion is larger. So these states 
would have an energy some MeV larger than respectively 
2.957 GeV and 6.370 Gev. 

A more quantitative computation of the masses, sizes, 
and decay rates of the proposed heptaquarks, includ- 
ing coupled channels and exact three body computations 
will be done elsewhere. I expect that the most relevant 
contributions that remain to be included in this frame- 
work are the solution of the full three body relativistic 
Fadeev equations, and the inclusion of the medium range 



interaction. The medium range interaction, which in nu- 
clear physics is described by the sigma meson exchange, 
is equivalent to the coupling to channels with multiple 
pions. 
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